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Abstract. On the verge of the centenary of dimensional analysis (DA), we present 
a new matrix generalisation of the Buckingham Theorem. The proof is based on 
a solution we have found for inverting non-square block matrices and gives rise 
naturally to a new pair of transforms - the similarity transform (S) that converts 
physical dimensional data into dimcnsionless space and its inverse (S ). Although it 
is well known that DA: a) reduces the number of free parameters, b) guarantees scale 
invariance through dimensional homogeneity and c) extracts functional information 
encoded in the dimcnsionless grouping of variables, scientists seem to be unaware 
that the scaling laws provided by DA are degenerate and therefore not unique. We 
demonstrate that the inverse transform S is responsible for the non-uniqueness and 
we show how reference to observational data is sufficient to to break the degeneracy 
inherent in transforming back to dimensional (physical) space. As an example, we 
show how the underlying functional form of the Planck Radiation Law can be deduced 
in only a few lines using the matrix method and without appealing to first principles 
- thus demonstrating the possibility of a priori knowledge discovery. It is hoped that 
the proof presented here will provide new impetus to the pursuit of inverse problems 
in physics. 
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1. Introduction 

The outstanding advances in experimental and theoretical physics in the second half of 
the 20 th Century have meant that we have been able to refine and improve the precision 
of the laws of physics (as evidenced by for example 15 decimal place accuracies for some 
atomic constants. However, the complexity of physical models has grown drastically at 
a rate which is proportional to the rate of accumulation of high resolution observational 
data meaning that exact, asymptotic and, in many cases, even approximate analytic 
expressions for the greater number of physical effects we are now able to include in our 
models, are difficult to deduce. As a result, many of us are resorting to systems-based 
modeling using for example holistic approaches such as complexity theory, or else are 
simulating interactions numerically using cellular automata, finite element analysis or 
Monte Carlo methods in an attempt to gain qualitative insight into complex natural 
phenomena. Furthermore, despite the giant technological leaps, observational errors, 
particularly in the field of astrophysics are of the order of 10%, the empirical relations 
we are obtaining are still very approximate (order of magnitude) and in their infancy 
|32j . In budding areas of research, DA is often used to obtain functional dependencies 
between the physical variables and the resultant scaling relations are then used to esti- 
mate the value of a given parameter. Although often, "back-of-the-envelope" in nature, 
the approach helps to constrain the problem. In this paper, we generalize the Bucking- 
ham Theorem upon which DA is based and uncover an important problem at its crux - 
that of the non-uniqueness of scaling laws. Though this may appear to be catastrophic, 
we present a methodology for its resolution. 

2. Dimensional analysis (DA) 

DA appeared early in the thoughts of physicists [9J, and since then has been studied 
and used by some of the most famous of them (for example, Maxwell [19], Rayleigh [26] 
and Einstein [8]). The first mathematical formulation of DA was published exactly a 
century ago [22] but it is often associated with the first general exposition of the ideas 
published in Nature by Buckingham [I] (see [25] for a discussion of its origins) . As such 
the fundamental theorem is often referred to as the Buckingham Theorem which has 
since been further developed in many textbooks [10], [15], [29] . 

Although scaling relations obtained by DA are prevalent in physics and astrophysics, 
a citation search in NASA ADS of all physics, astrophysics and electronically archived 
papers, for example, shows that only 35 refereed articles have ever cited Buckingham's 
original paper and all have failed to address the problem of the uniqueness of dimen- 
sionless numbers. This means that scaling relations used by or derived in these papers 
require great caution. The problem of uniqueness has been partly addressed in the 
pioneering work of Stephan Rudolph and co-workers [13], [27], [28] from the engineering 
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design community. In this paper, we show using a matrix generalisation of DA precisely 
where the degeracy comes from and how it may be addressed. Before doing so, we first 
clarify some of the lesser-known pitfalls associated with dimensionless numbers that, as 
we shall see, are the key to an efficient pre-processing of data. 

2.1. Dimensionless numbers 

A dimensionless number (or more precisely, a number with the dimensions of "1") has 
no physical units. It is typically defined as the ratio of quantities having units of iden- 
tical dimension whose units then cancel. Whenever one measures a physical quantity 
it is with reference to a dimensional standard so that ultimately, we always work with 
dimensionless numbers and ratio scale, physical quantities. It is essential that the units 
are the same in both the numerator and denominator, such as m/m to avoid errors 
associated with unit conversions (e.g. cm/m and mm/cm are also dimensionless but 
introduce scaling factors of 1/100 and 1/10 respectively with respect to m/m). Already, 
here, we can see a degeneracy creeping in. We can, in fact, construct an infinite set of 
such length ratios that would all be dimensionless but which would differ by a constant 
of multiplication. Furthermore, a physical quantity that may be dimensionless in one 
system of units, may not be dimensionless in another system of units. For example, in 
the non-rationalized CGS system of units, the unit of electric charge (the Stratcoulomb) 
is defined in such a way so that the permittivity of free space eo = 1/4-7T whereas in the 
rationalized SI system, eo = 8.85419 x 10~ 12 F/m. In systems of natural units (e.g. 
Planck units or atomic units), the physical units are defined in such a way that several 
fundamental constants are made dimensionless and set to " 1" thus removing these scal- 
ing factors from equations. While this is convenient in some contexts, abolishing all or 
most units and dimensions makes practical physical calculations more error prone PQ, 
especially, as we will see, when the goal is to pre-process data and extract empirical 
relations. The next section describes Buckingham's simple and elegant theorem con- 
cerning dimensionless numbers. 

2.2. The Buckingham ix -Theorem 

In the physical sciences, the two sides of any equation must be commensurable or have 
the same dimensions (the principle of dimensional homogeneity) for all valid functional 
relationships [5j. This led Buckingham to prove [I] that, from the existence of a complete 
and dimensionally homogeneous function, / of exactly n physical variables {x n } € 9ft + , 



there exists a corresponding functional relationship F of m < n dimensionless numbers 




(1) 




(2) 
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where m = n — r is reduced by the number (r) of dimensionally independent variables in 
{x n }. The restriction to positive values of the variables can be satisfied universally by 
suitable coordinate transformations. The Buckingham Theorem holds for all dimension- 
ally correct equations of physics [23], [10] . The 7r-groups are dimensionless numbers or 
"similarity numbers" which are the scaling invariants of a physical problem. Although 
simple, in essence, Buckingham's Theorem is extremely profound. From our observa- 
tions of reality we categorise data into homogeneous sets that we identify with different 
physical variables {x n }. By constructing a qualitative model relating the variables, we 
find p different various physical properties f(xi, ■ ■ • , x n ) p . In the context of a chosen sys- 
tem of units (e.g. SI) and frame of reference (e.g. coordinate system or inertial frame), 
we are able to offer a quantitative description of reality f(xi, ■ ■ ■ ,x n ) = 0. Finally, in 
order to correctly evaluate reality on all scales, we must find the mimimal and invariant 
description in terms of dimensionless numbers F(tti, ■ ■ • , 7r m ) = 0. It is at this point that 
we see the importance of DA - only once the function F is found have we successfully 
decoded the empirical law describing the data. 

2.3. A Matrix Generalization of DA 

DA provides a procedure to generate dimensionless numbers from a list of relevant 
physical variables together with their respective dimensions. Several authors have 
heuristically arrived at matrix expressions to calculate dimensionless numbers |29j.[27] 
but a full mathematical proof has, until now, been lacking. The reason appears to be 
two-fold: 

(i) for simple problems, a dimensionless combination of variables can often be found 
by inspection (e.g. a ratio of lengths), 

(ii) a generalised matrix inverse for the resultant non-square, degenerate matrices 
involved (see below) has, until now, only been calculable numerically and therefore 
approximately. 

It is possible that the failure of many authors to address the problem of the non- 
uniqueness of dimensionless numbers may be traced back to the absence of such a proof, 
which we present below. In what follows, we use the notation employed in |29|. 

The dimensional equations of the physical variables in equation 1 (with " [ ] " repre- 
senting "dimensions of) can be expressed as: 



where the {dj,j = l,...,r} are the base dimensions having exponents a^. In the 
SI system of units, there are 7 distinct base dimensions: mass M[Kg], length L[m], 



N 



d^ l d^\..d^ 
d« 21 d% 22 ...d« 2r 




d" nl d% n2 ...d* 
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time T[s], temperature 9[K], electrical current I [A], concentration N[mol] and light 
intensity J[cd]. So, for example, a pressure P which is measured in N/m 2 = [kg 1 m~ 1 s~ 2 ] 
would have a dimensional equation of the form [P] = M l L~ l T~ 2 0° 'I°mol J° . In this 
case, the exponents of the dimensions of pressure are given by the row vector 
(1, —1, —2, 0, 0, 0, 0). The matrix formed with the exponents «y of the dimensions dj of 
all physical variables Xi relevant to a problem, we call the dimension matrix D: 



D 



an Oii2 . . . ai r 

a 2 l a 2 2 • • • «2r 



(3) 



D has n rows for the number of physical variables and r columns for the distinct base 
dimensions. To satisfy dimensional homogeneity, we need to find the exponents e« that 
solve the equation, 

[ x i' x 2 2 --- x 'n\ = d?2 . . .d q r T , 
with {qj, j = 1, •••,/*} being the sought combinations of fundamental dimensions (to 
be set equal to zero for dimensionless numbers). The condition, q 1 — q 2 • • • q 3 ■, — then 
leads directly to dimensional homogeneity: 

K*2 2 -< n } = 1- (4) 
A dimensionless (dimension " 1" ) n group is therefore given by, 

Tr = x\ 1 X? (5) 

Our task here then is to find the exponents €j. We proceed as follows. The dimensional 
equations for the physical variables [xj], when raised to the power 6j become, 

[ x ^] = (d^d^.-.d^y 1 
[ x %\ = (d^d^ 2 ...d^y 2 

[ x ^\ = (d% nl d% n2 ...d? nr y n . 
Evaluating the product 117=1 l x T] an( i grouping the base dimensions dj we obtain, 

Irr^lrr^ rr^n} / JOniei+a2ie2 + ---+a n l€ n 

l x l x 2 ■■~ L n \ ~ "l 

X ^«12«l+Q!22«2 +— +an2«n 
X Q r ^-^^\+OL2r<^2 + ---+OL rir <L ri 

Associating qj with the exponents of the base dimensions dj then we arrive at the 
following linear system of equations, 

anei + a 2 ie 2 + . . . + a n ie n = qi 
ai 2 ei + a 22 e 2 + . . . + a n2 e n = q 2 



ai r €i + a 2r e 2 + . . . + a nr e n = q r , 
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which in matrix form are given by, 

D T e = q. (6) 
The €j can then be obtained from the inverse dimension matrix, 



D 



T 



provided that 
a single dimension 



P 1 ') q. (7) 

7^ 0. The vector e containing the ej then allows identification of 
ess number. However, we know from Buckingham's Theorem that 
there are exactly m such groups. We therefore require a solution vector for each of the 
7r-groups. Thus, in equation (GO), the vectors e and q must become matrices e and q of 
dimension (n x m) and (r x m) respectively. The result is that dimensional homogeneity 
will then lead to the set of (j — 1, • • • , m) dimensionless groups iij given by, 

71 j = x\ 1J X2 J ■ ■ ■ X e n "i . (8) 

Our next task is to invert the dimension matrix D T . However, it is non-square and 
degenerate since the number of physical variables is always greater than the number of 
related base dimensions |5]. The generalised inverse of a non-square matrix is far from 
trivial to find and, until now, only numerical approximations have been available [21J. 
However, for a special class of inverses |7] we know that the matrix equation Ax = b for 
a non-square (n x m) matrix A, has a solution iff AA + b = b where the (m x n) matrix 
A + is the Moore-Penrose inverse of A such that AA + A = A. In a second remarkable 
paper published in the same year [16J , it was shown that the Moore-Penrose inverse A + 
is rarely equal to the generalised matrix inverse A" 1 . The matrix solution provided by 
[29] obtains the Moore-Penrose inverse. Although A + can always be found numerically 
by SVU decomposition for well-defined problems, our interest here is the closed (exact) 
algebraic form. 

Rather than seeking the Moore-Penrose inverse and then attempting to show its gen- 
erality for the infinite set of solutions, we have been able to find a closed form for the 
generalized solution directly. The proof is short and proceeds as follows. We know from 
Buckingham's Theorem and the principle of dimensional homogeneity that each dimen- 
sionless 7r-group is constructed from an independent physical variable multiplied by a 
suitable combination of the dependent physical variables [I]. We start then, by parti- 
tioning the dimension matrix D T into two inner block matrices: A T for the dependent 
variables and B T for the independent variables, 



D 



T 



B T A 1 



(9) 

The rank r of the matrix A T is equal to the number of base dimensions dj involved in 
the problem. It is known [21] that if X = A -1 is any matrix satisfying, 

AX A = A, (10) 

then the linear system Ax = b has a solution iff, AXb = b. Such as system has a general 
solution, 

x = Xb+(l-XA)y, (11) 



100 Years of dimensional analysis: new steps toward empirical law deduction 



7 



with y arbitrary. In order to calculate the complete set of it groups, we need to solve 
the matrix system equivalent of equation (jSJ), 

D T i = q 



such that, 



D 



T 



-1 



(12) 
(13) 



Here, e is a matrix of size {n x m), q is a matrix of size (r x m) and D T is a block matrix 
of size (r x n) partitioned into the (r x m) matrix B T and the (r x r) invertible matrix 
A T . Taking into account that from equation (flOj) . 







B T , A T 



and that, 



5 T A T 



A T ~ 



q = iq = q, 



then, ( fill) provides the general solution for equation ( fl3l ): 





A T ~\ 



q+ [I 



+ 



A?" 1 



I 




/ 





-i T 1 B T 




2/1 

y~2 



2/1 



A T g — A T B T y~i 
q - B T y~\ 



AT" 1 



(14) 



Here, y~i is an arbitrary matrix of size (m x m). Equation (|T4"|) gives the whole set 
of solutions of the matrix system ( Tl2l) . Here, we are concerned with obtaining the 
dimensionless numbers irj for which we need to set q = 0. Inserting this into equation 
(|T4|) we see that, in order to satisfy equation (Tl2l . the arbitrary matrix yj must be equal 
to the identity matrix J since only then do we find that 



D T ~e 



B T A 1 



B T I — 
B T I — 



A T (A T ) 
IB T = a 



-A T X B T 

-1 ~rr 

B T 
= 0. 
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Figure 1. Construction of the dimensionless matrix from the dimensional one using 
rank-preseving operations 



The condition q = is a special case for which the arbitrary matrix y~\ becomes specfic 
rather than arbitrary. This is exactly the case we are seeking to determine dimensionless 
combinations of the physical variables. The solution to our problem then, is given by, 

' (15) 



e(q = 0) 



-A T 



1 B T 



The solution matrix e has n rows, one for each of the physical variables X{ and m columns, 
one for each of the dimensionless numbers ttj. The identity matrix / in the upper half of 
e guarantees that only one independent variable at a time is combined with dependent 
variables in constructing the 7t-groups. Furthermore, as we might expect, the solution 
matrix e is a block matrix having an upper-diagonal form. In the case of the inversion by 
Guassian- Jordan elimination of a non-singular square matrix A , then problem matrix 
is augmented with the identity matrix and inversion proceeds as follows: 



A 



A s I / 



I I (A s ) 



-i 



For our non-square partitioned matrix D T then, in obtaining the inverse matrix, only 
the symmetric part of B is effectively augmented by the identity matrix; the semantic 
structure of Guassian Elimination is maintained. Figure [1] shows the overall process. 
The dimension matrix D, (left), is formed from the exponents oty of the base dimensions 
dj of the relevant physical variables Xi\ one row for each variable and one column for 
each base dimension. The rank r of the dimension matrix is equal to the number of 
base dimensions involved and is used to partition D (centre) with the block matrix B 
containing the exponents of the independent variables and the block matrix A contain- 
ing the exponents of the dependent variables. The solution matrix e is then obtained 
and comprises an upper square identity matrix / and a lower sub- matrix — (A T )~ l B T . 



The various matrices used in the proof include the square matrices I mm , A rr and the 
non-square matrices: D nr , DT nr , (D T )~^, e nm , q rm , B mr and (A T )~ 1 Bj m . The math- 
ematical complexity of the problem of evaluating the dimensionless 7r-groups therefore 
lies in the calculation of the inverse matrix (A T )~ 1 of dimension equal to the rank r of 
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D. The condition for (A T )~ 1 to exist is that the determinant A T ^ 0. This mathemat- 
ical condition corresponds to the physical condition that we have correctly identified 
in x n the complete list of relevant physical variables required to satisfy dimensional 
homogeneity. Had we missed out a crucial variable, then the determinant would equal 
zero producing a singular and undefined inverse. We thus, have an additional check on 
the correct construction of the problem in dimensional space. Since the dimensionless 
groups form a reduced set (there being m = n — r of them), they minimise the number 
of free parameters. The 7r-groups then provide the most compact formulation possible 
for the evaluation of physical laws. In the next section we show how the matrix method 
gives rise naturally to a pair of transforms: the Similarity Transform S and its inverse 
S~ x that allow us to reveal the origin of the problem of non-uniqueness associated with 
scaling laws. 



2.4- Similarity Theory and the Uniqueness Problem 

We have described how functions of the type / (x\, ■ ■ ■ ,x n ) = or equivalently, 
x n = / (xi, • • • i Xn-i), serve as quantitative models. The philosophical principle of 
causality [6] is guaranteed when, for a complete list of physical variables x n , the presence 
of cause (xi, ■ ■ ■ , x n -i) always results in effect x n . Any question about the behavior of 
the physical object can be answered by the rule /, otherwise / (xi, ■ • • , x n ) = would not 
represent the complete and correct physical model. This means that F (71*1, ■ • ■ , 7r m ) = 
represents the similarity transform of / under which the physical content of the rule 
F remains invariant [28]. The similarity transform S : X 1— ► II from the space X of 
physical variables {x{\ to the dimensionless space II of dimensionless variables {7^} from 
equation (jSJ) is given by, 

n 

II (16) 

i=i 

with j G [l,77i]. S represents the transform of a space 9ft n onto a space 5i m where 
m = n — r corresponds to a dimensionality reduction leading to different points in X 
being mapped to the same point in II. From the Buckingham Theorem for the functional 
relationship between 7r groups, 

F(7Ti, ■ ■ ■ ,7T m ) = 0, 

any single 7T group can. analogously, be expressed as a causal function F of the remaining 
groups: 

ttj = F(Tr^) (17) 

with k < m. The 7r groups are calculated from the physical variables x n using the 
Similarity Transform S (equation [T6i) . Since each ir group contains one and only one 
independent variable, we can find the inverse transform using their negative exponents 
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in the product functio such that the other independent variables cancel out in the j th 
group, 

r 

e (i+m)j 

1=1 

Then, substituing for ttj from equation (1171) . the inverse similarity transform S" _1 is given 
by, 

r / n \ 

S- 1 : Xj = J] ^ mh F U k = J] ^ • (18) 

1=1 \ 1=1 / kj^j 

This then gives the scaling laws for the independent variables Xj expressed as a function 
F of the dependent variables combined in 7r groups. 

The inverse similarity transform 5 1-1 : II i— > X of the dimensionless space II to the 
dimensional space of physical variables X is not unique since it maps a space 3? m onto 
a larger space 3? n ; there being only m dependent variables. A single point p in X- 
space with (x\, ■ ■ ■ ,x n ) is mapped via the S transform onto its corresponding point 
(•7Ti, • • • , Tc m ) p in dimensionless space II. The inverse transform S" -1 leads to the whole 
set of all completely similar points (xi, ■ • ■ , x n ) q=1 ... ^ in X defined by the similarity 
condition [tti, ■ ■ ■ , ^m) v=consta nf ^ n e function F in Il-space is simpler than its coun- 
terpart / in dimensional space X; being of lower dimensionality. Mathematically, the 
7Tj=constant define a hypersurface in X that contain all completely similar points and 
reflect the fact that there exists an infinite set of dimensionless numbers for a given 
physical problem [28]. This can be seen by taking, for example, the dimensionless ratio 
7i l = h/h of two lengths l\ and li. When l\ = 3m and 1% = or when l\ = 6m and 
I2 = 8m, then in both cases ttl = 0.75. It is the non-uniquness of the inverse mapping 
II 1— > X that is responsible for the non-uniqueness of F and hence of scaling laws de- 
duced from the inverse transform S" _1 . Figure [2] shows a visual representation of how 
the asymmetry between / and F results in non-uniqueness of the dimensionless groups. 
It is this non-uniqueness that propagates back to scaling laws and the identification of 
empirical laws. 

DA provides the mechanism for mapping a problem into a space where dimensional 
reduction can be exploited. Proper identification of the function F then allows for a 
full problem evaluation obtained by inverse transformation back to dimensional space. 
Fourier transforms, Laplace transforms and the like all take advantage of the algebraic 
simplicity in the transformed space before inverse-transforming back to the problem 
space. However, it is the ambiguity that arises from the inverse transform that is ul- 
timately responsible for the problem of non-uniqueness when the dependent variables 
are expressed in terms of the tt groups. As an illustration, in the next section we apply 
DA to the problem of Planck Radiation to show a) the ease of application of our ma- 
trix method, b) how the problem of degeneracy manifests itself, c) how, despite this a 
priori knwoledge discovery is still possible, and d) how reference to data can break the 
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Figure 2. Similarity transforms: dimensionality expansion 7r 1 of case p onto all 
completely similar cases q = 1, • • • , oo: adapted from [28] 



degeneracy and provide unique scaling laws. 
3. The Planck Radiation Law 

For Planck radiation, we know that the spectral intensity u is related to the 
electromagnetic radiation emitted at different frequencies v = c/X from a black body at 
a temperature T. Since the radiation is generated by oscillators having energy k B T per 
degree of freedom and is carried by photons travelling at the speed of light c in packets 
of energy hu, we expect that the total integrated spectral intensity will be a function of 
the following list of physical variables: 

f(u,X,h,c, k B} T) = (19) 

The dimensional equations of the physical variables are: 



[u] 


= M 1 L- 2 T- 2 d® = ML~ 2 T~ 2 


[A] 


= M°L 1 T°6° = L 


[h] 


= M 1 L 2 T~ 1 9° = L = ML 2 T 


[c] 


= M°L 1 T~ 1 9° = IT' 1 


[k B ] 


= m 1 l 2 t- 2 6- 1 


[T] 


= M^L^O 1 = 0, 
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and give rise to the dimension matrix, 



D 





M 

1V1 


T 

±J 


T 


a 
u 


u 


1 


-2 


-2 





X 





1 








h 


1 


2 


-1 





c 





1 


-1 





k B 


1 


2 


-2 


-1 


T 











1 



(21) 



The next step involves seperation of the elements of D T into two block matrices: A T for 
the dependent variables and B T for the independent variables. Since there are n = 6 
physical variables in equation ( Tl9l) and r = 4 base dimensions, the number of 7r groups 
will be m = n — r = 2. There are therefore r = 4 dependent variables and m = 2 
independent variables. We are looking for an equation for the spectral intensity u and 
so this will be one of the two independent variables; the other is the wavelength A. As a 
check on the completeness of the list of physical variables, we calculate the determinant 
\A T \ = 1 ^ as required for dimensional homogeneity and for the invertibility of A. We 
can now partition the dimension matrix D and take the transpose, 



D 1 



B T A T 







1 

2 



1 
1 

-1 -1 











1 



(22) 



The matrix product — (A ) 1 B T then gives the solution matrix, 



'B 1 





7Tl 


vr 2 


u 


1 





A 





1 


h 


4 


-1 


c 


4 


-1 


k B 


-5 


1 


T 


-5 


1 



(23) 



The m 



2 different 7r-groups are then obtained from the similarity transform S ( 116 

w/i 4 c 4 



711 



712 



u 1 X°h i c i k- B b T- 



'c-'kiT 1 



k%T 5 
XkT 



^ he 

The inverse similarity transform S 1 fj 18)) with x\ 

hcF(n 2 



u 



k%T 5 



F(ti 2 ) 



7T 5 




(24) 



(25) 



/i 4 c 4 v ' A 5 

Comparing the result of DA with the exact expression for Planck's Law for the spectral 
energy density per unit wavelength integrated over all solid angles: 
_ he 1 



u 



•57T 



e fhc/Xk B T\ _ I ' 



(26) 
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we see that remarkable progress has been made using a priori knowledge discovery alone. 
However, DA cannot provide the values of constant multipliers such as the "87r" or the 
actual form of the function F = l/i^e 712 — 1). And what about non-uniqueness? Since, 
7T2 is dimensionless, then so too are any algebraic combinations of 7r 2 such as tt^ , \ or 
1 + 7T2 and so forth ad infinitum. The question we are forced to ask now is the following: 
faced with this problem of uniqueness and the degeneracy created in the function F, 
how are we to find objective and unique dimensionless scaling relations? In particular, 
empirical laws published in the literature that do not resolve the problem of degeneracy 
should be treated with great caution. However, the information encoded into actual ob- 
servational data is both sufficient and necessary to uniquely and unambiguously deduce 
F and break the degeneracy. In a recent paper [32] we applied an evolutionary genetic 
network to a sample of 98 young galaxies where we obtained the lowest error empirical 
relation for their chemical abundances to date. Given a real data set, there will be one 
and only one best fit function F which can be used to break the degeneracy inherent in 
the purely deductive method of DA presented in this paper. 

4. Conclusions 

We have generalized DA using new matrix theory that is easy to apply. We have 
shown how the determinant check ensures a complete specification of the problem and 
that the theory gives rise to a similarity transform S and its inverse S' 1 that depends 
on a non-unique and degenerate functional relationship F between the dimensionless tt 
groups. The example of Planck's Law shows that DA is capable of substantial functional 
knowledge discovery that can be used as an intermediate step prior to data analysis, 
but that a full and unique solution is only possible with a correct identification of F by 
using observational data to break its degeneracy. We hope that the simplicity of the 
matrix method and its ease of application will help pave the way for a new approach to 
the search for empirical laws and scaling relations from data. 
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